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LONG-TERM  GOALS 

Terrestrial  runoff  and  river  input  dominates  urban  pollutant  loading  rates  Often  draining  directly 
into  the  surfzone,  this  pollution  degrades  surfzone  water  quality,  leading  to  beach  closures 
(e.g.,  Boehm  et  al.,  2002),  increases  health  risks  (e.g.,  diarrhea  and  upper  respiratory  illness) 
(Haile  et  al.,  1999),  and  contains  both  human  viruses  (Jiang  and  Chu,  2004)  and  elevated  levels 
of  fecal  indicator  bacteria  (Reeves  et  al.,  2004).  Surfzone  mixing  processes  disperse  and  dilute 
such  (and  other  types  of)  pollution.  On  smaller  length-scales  (smaller  than  the  water  depth), 
breaking-waves  and  bed-generated  turbulence  mix  tracer.  However,  field  surfzone  observations  of 
turbulence  previously  have  been  extremely  scarce,  and  much  about  surfzone  small-scale 
turbulence  is  not  known.  On  larger  scales  (10-100  m),  horizontal  dispersion  is  driven  by  surfzone 
eddies  and  meanders  associated  with  shear  waves  (Oltman-Shay  et  al.,  1989)  or  finite  breaking 
crest  length  (Peregrine,  1998).  Understanding  the  small  and  large  length-scale  mixing  processes 
important  to  predicting  the  fate  (transport,  dispersal,  and  dilution)  of  surfzone  tracers  whether 
pollution,  bacteria,  larvae,  or  nutrients. 

OBJECTIVES 

The  scientific  objective  is  to  improve  understanding  and  modeling  of  dispersion  of  tracers 
(pollution,  fecal  indicator  bacteria,  fine  sediments)  within  the  nearshore  (a  few  100  m  of  the 
shoreline)  and  especially  within  the  surfzone  where  breaking  waves  intensify  mixing  processes 
and  drive  strong  mean  currents.  Here  the  focus  is  on  two  components  of  the  analysis  of  the  HB06 
experiment  performed  in  Fall  of  2006.  The  first  is  analysis  of  surfzone  cross-shore  dye  tracer 
dispersion  (Clark  et  al.,  2009b).  The  second  is  studying  the  small-scale  turbulence  in  the  surfzone 
due  to  breaking  waves.  Other  HB06  efforts  include  analysis  of  surfzone  drifter  dispersion 
(Spy dell  et  al.,  2009)  and  phytoplankton  patchiness  (Omand  et  al.,  2009b). 

APPROACH 

HB06  Dye  Dispersion 
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Previous  surfzone  dye  observations  have  been  limited  to  discrete  shoreline  or  single  swimmer 
bottle  samples  without  measuring  the  waves  and  currents  ( e.g Inman  et  al.,  1971;  Clarke  et  al., 
2007)  A  newly  developed  ONR-funded  GPS  equipped  jetski  based  sampling  system  ( Clark  et  al., 
2009a)  was  used  to  create  near-synoptic  maps  of  surfzone  dye  (Rhodamine  WT)  on  six  days 
(denoted  R1-R6)  during  HB06  ( Clark  et  al.,  2009b).  Dye  tracer,  released  relatively  close  to  the 
shoreline,  was  transported  alongshore  by  the  alongshore  current  and  dispersed  in  the  cross-shore 
in  a  manner  resembling  a  wall-bounded  turbulent  plume,  analogous  to  a  smokestack  plume,  with 
axis  parallel  to  the  shoreline  ( Clark  et  al.,  2009b). 

Surfzone  Turbulent  Dissipation  Rate: 

The  vertical  structure  of  turbulence  in  the  surfzone  is  of  interest.  Both  breaking  waves  and 
near-sea-bed  shear  are  possible  sources  of  turbulence.  Here  a  key  turbulence  statistic,  the 
turbulent  dissipation  rate  e  is  estimated  from  Acoustic  Doppler  Velocimeters  observations 
following  Feddersen  et  al.  (2007). 

WORK  COMPLETED 

•  Spydell  et  al.  (2009)  was  submitted  to  JGR  Oceans  in  Jan  2009  and  was  published  this 
summer.  This  manuscript  reported  on  HB06  drifter  observations 

•  Omand  et  al.  (2009)  and  Clark  et  al.  (2009a)  have  both  been  published.  These  manuscripts 
dealt  with  the  technical  details  of  making  Chlorophyll  and  Rhodamine  WT  dye 
observations  in  the  surfzone. 

•  A  manuscript  ( Clark  et  al.,  2009b)  of  the  HB06  dye  dispersion  studies  has  been  submitted 
to  JGR  Oceans. 

•  A  manuscript  ( Omand  et  al.,  2009b),  reporting  on  the  evolution  and  dynamics  of  a 
nearshore  red  tide  observed  during  HB06,  is  near  submission. 

•  A  manuscript  intended  for  J.  Atmospheric  and  Oceanic  Tech.  (Feddersen)  is  almost  ready 
for  submission.  This  manuscript  deals  with  the  methods  of  analyzing  Acoustic  Doppler 
Velocimeter  data  for  estimating  the  turbulent  dissipation  rate,  manuscript  looks  at  surfzone 
drifter  dispersion  from  Huntington  Beach. 

•  Preparations  for  the  IB09  (Imperial  Beach  CA  in  Sept-Oct  2009)  experiment  have  been 
completed  and  the  IB09  experiment  is  currently  under  way. 


RESULTS 

HB06  Experiment 

Observations  were  collected  from  15  September  to  17  October  2006  (800  hours)  at  Huntington 
Beach  CA  ,a  site  with  chronic  water  quality  problems.  A  cross-shore  transect  of  co-located 
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Distance  from  the  Mean  Shoreline  (m) 

Figure  1:  HB06  cross-shore  instrument  transect.  The  circled  numbers  indicate  instrumented  frames 
with  acoustic  Doppler  Velocimeters  (ADV)  and  pressure  sensors.  The  horizontal  red  dashed  lines 
indicate  the  typical  tide  range. 


pressure  sensors  and  acoustic  Doppler  Velocimeters  was  deployed  spanning  160  m  out  to  4  m 
mean  water  depth  (Fig.  1).  The  tide  range  was  nominally  ±1  m.  The  data  was  sampled  at  8  Hz. 
The  ADVs  sampled  between  0.5- 1.0  m  above  the  bed.  The  cross-  and  alongshore  coordinate  are  x 
and  y,  respectively.  The  mean  water  depth  is  given  by  h  and  the  vertical  coordinate  is  z  with 
z  =  0  m  at  the  sea-bed.  The  distance  below  the  mean  sea  surface  is  z!  —  h  —  z.  At  each  of  the 
frames,  hourly  estimates  of  significant  wave  height  //sig,  cross-shore  energy  flux  F  =  Ecg,  mean 
alongshore  current  vm,  and  turbulent  dissipation  rate  e  were  estimated. 


HB06  Dye  Dispersion  Obsen’ations 


The  mean  cross-shore  tracer  profiles  D(x,  yj )  average  over  stirring  and  meandering  (Fig  2)  On  all 
releases  except  R5,  the  initially  narrow  D(x,  yj)  profiles  disperse  across  the  surfzone  and  peak 
concentrations  decrease  with  downstream  distance  (increasing  y)  from  the  source  (Fig  2).  During 
most  releases  (R2,  R3,  R4,  R6),  mean  profiles  D(x,  yj)  are  shoreline  attached  with  maxima  at  or 
near  the  shoreline  (Fig  2).  An  exception  is  release  R1  which  has  maxima  in  the  mid  to  outer 
surfzone  (Fig  2a),  which  is  likely  the  result  of  dye  released  in  the  mid-  to  outer-surfzone 
(iCo  =  —55  m),  and  the  shortest  yr  Release  R5  was  sampled  far  downstream  of  the  the  dye  source 
(y  >  240  m)  where  tracer  had  already  spread  across  the  surfzone. 

A  simple  model  for  a  shoreline  attached  tracer  plume  is  ( Clark  et  al,  2009b) 


D(x,tp) 


Qo 


\J  47TKXa;fj 


exp 


—  lx 


—  2  \ 


4/CxX^P 


(1) 


where  nxx  is  the  cross-shore  diffusivity,  tp  =  y/V  is  pseudo-time  (V  is  the  surfzone-averaged 
mean  alongshore  current),  and  Q0  is  the  normalized  flux  of  (known)  tracer  input  at  the  source. 
This  model  (1)  has  the  cross-shore  tracer  half- width  <7gUrf  =  nxxtp  growing  linearly  with  time.  For 
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Figure  2:  Mean  cross-shore  tracer  concentration  D(x,  y-j)  versus  x  for  releases  (a)  Rl,  (b)  R2,  (c) 
R3,  (d)  R4,  (e)  R5,  and  (f)  R6.  Colors  indicate  different  downstream  alongshore  distances  y  (see 
legends  in  each  panel).  The  vertical  dashed-gray  line  represents  the  outer-limit  of  the  surfzone. 


5  of  the  6  releases  (R5  kxx  could  not  be  estimated),  surfzone  kxx  (together  with  error  bars)  is 
estimated  from  the  slope  of  the  least  squares  cr2urf  versus  tp  fits  (e.g.,  for  R6,  see  Fig.  3).  For  Rl, 
the  Kxx  estimation  method  is  modified  to  account  for  non-shoreline  attached  conditions.  Only 
transects  where  tracer  is  considered  surfzone-contained  are  used  in  the  fits  (e.g.,  solid  black 
symbols  in  Fig.  3).  Estimated  kxx  range  from  0.5-2. 5  ±  0.62  m2  s-1.  For  R6,  once  the  tracer 
saturates  the  surfzone,  the  plume  width  stops  growing  (e.g.,  open  symbols  in  Fig.  3).  However, 
this  was  not  generally  observed  as  jetski  transects  were  not  typically  driven  at  such  large  y  (or 
large  tp).  Proposed  dye  release  experiments  will  include  such  large  downstream  sampling. 

The  simple  model  (1)  predicts  for  maximum  tracer  concentration  -D^ax  =  Q0/ nxxtv)] for 
(near- )  shoreline  releases.  The  estimated  surfzone  kxx  and  the  known  injected  dye  flux  (Q 0)  are 
used  to  predict  D ^ax  for  shoreline- attached  releases  (R2,  R3,  R4,  and  R6)  when  tracer  was 
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Figure  3:  ofurf  (±  error-bars)  versus  tp  for  release  six  (R6).  Black  symbols  indicate  observations 
when  tracer  was  surfzone-contained.  Fits  to  estimate  kxx  are  indicated  with  the  dashed  gray  line. 


surfzone-contained.  The  observed  Dmax  are  consistent  with  predicted  P^.AY  (Fig.  4)  with  high 
skill  (0.72),  although  the  R3  and  R6  have  large  errors  at  the  first  downstream  transect.  The 
linear  growth  of  ofurf  and  the  skill  in  predicting  D^ax  indicates  that  the  simple  model  (1)  well 
represents  the  surfzone  mean  tracer  plume  structure  as  it  is  advected  downstream.  The  initial 
stages  of  dye  dispersion  (the  near  field)  and  the  far  downstream  are  poorly  understood. 

Two  kxx  parameterizations  for  tracer  dispersion  were  examined  by  Clark  et  al.  (2009b).  The  first 
is  a  breaking-wave  (bore)  induced  mixing  scaling  where  kxx  ~  ( Feddersen ,  2007; 

Henderson ,  2007)  ( Hs  is  wave  height  and  Tm  is  wave  period).  In  drifter  model  simulations 
( Spydell  et  al.,  2009),  horizontal  rotational  velocities  (i.e.,  vortical  flow)  generated  by  finite  crest 
length  breaking  ( Peregrine ,  1998)  or  shear  instabilities  of  the  alongshore  current 
( e.g .,  Oltman-Shay  et  al.,  1989)  were  found  to  be  a  primary  mixing  mechanism.  A  mixing-length 
scaling  (e.g.,  Tennekes  and  Lumley,  1972)  was  examined  using  a  surfzone  width  Lx  length-scale 
and  infragravity  horizontal  rotational  Vrot  velocities  for  a  velocity  scale,  i.e., 

K'xx  (2) 

where  a  is  a  constant  expected  to  be  0(1)  but  <  1.  The  incident  Hs  and  Tm  are  used  to  test  the 
bore  induced  kxx  scaling.  The  observed  surfzone  Vm1  (x)  is  estimated  following  Lippmann  et  al. 
(1999),  and  the  surfzone  averaged  Vm1  range  between  0.036-0.09  ms-1. 

Although  observed  kxx  generally  increase  with  IP  T~ 1  (Fig  5a),  the  skill  is  low  (r2  =  0.32)  and 
the  fit-slope  of  11.7  is  a  factor  6  larger  than  expected  for  bore-induced  dispersion  ( Clark  et  al., 
2009b),  suggesting  that  the  observed  cross- shore  dye  dispersion  is  not  dominated  by  the 
bore-mixing  mechanism.  The  surfzone  tracer  kxx  increase  with  Vrot  Lx  (Fig  3b)  and  the  linear 
best-fit  gives  an  r2  of  0.59,  slope  of  0.2,  and  near-zero  y-intercept.  The  high  r2  and  an  expected 
slope  <  1  (for  a  mixing-length  scaling)  indicate  that  rotational  velocities  (surf-zone  eddies)  play 
an  important  role  in  cross-shore  surfzone  tracer  mixing. 

However,  only  6  total  dye  releases  were  performed  at  a  single  beach  without  significant  variation 
in  the  wave  and  current  conditions.  In  addition,  the  nxx  with  overlapping  error  bars  (Fig  5)  are  not 
particularly  distinct  on  the  different  release  days.  Furthermore,  the  initial  short  time  (<  100  s  or 
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Figure  4:  Predicted  versus  observed  Dmax  for  shoreline  attached  releases  (R2,  R3,  R4,  R6)  and 
surfzone  contained  transects.  The  skill  is  0.72. 


the  near  field)  cross-shore  dispersion  of  the  tracer  was  not  consistently  observed  and  only  on  two 
release  days  were  the  long  time  (far  field,  y  >  500  m  downstream)  saturation  of  the  surfzone 
observed.  During  HB06,  dye  was  never  observed  more  than  700  m  downstream  of  the  dye  source. 
How  the  tracer  evolves  and  mixes  with  the  waters  offshore  of  the  surfzone  is  totally  not 
understood. 

Lastly,  the  jetski  only  measures  near-surface  dye  concentrations.  For  analysis  purposes,  the  dye  is 
assumed  to  be  well  mixed  (vertically  uniform)  within  the  surfzone.  Within  the  surfzone  where  the 
depths  are  shallow  (<  3  m)  this  assumption  does  not  appear  to  affect  cross-shore  integrated  dye 
moments  ( Clark  et  al.,  2009b).  However,  the  validity  of  this  assumption  is  unknown.  The  dye 
tracer  that  leaks  offshore  can  potentially  be  concentrated  at  the  surface  or  at  depth.  The  amount  of 
dye  that  leaks  offshore,  it’s  vertical  distribution,  and  the  mechanisms  that  drive  it  are  unknown. 

Surfzone  turbulence  studies 

The  turbulent  dissipation  rate  e  was  significantly  larger  (by  a  factor  of  10)  inside  the  surfzone 
(blue  curve  first  400  hrs)  versus  seaward  of  the  surfzone  (red  curve  in  Fig.  6),  indicating  the 
importance  of  wave  breaking  to  turbulence  in  the  surfzone.  The  surfzone  e  is  related  to  the 
incoming  Hsig  (lower  panel  Fig.  6)  with  larger  waves  leading  to  larger  e.  Furthermore,  there  is 
significant  tidal  modulation  of  surfzone  e.  At  lower  tides  when  the  ADV  is  closer  to  the  surface 
(and  the  source  of  breaking  wave  turbulence),  e  is  stronger.  In  addition,  at  later  times  (hours 
500-700),  frame  3  is  alternately  within  and  seaward  of  the  surfzone  as  the  tide  goes  up  and  down. 
The  frame  3  e  varies  strongly  as  wave  breaking  is  turned  on  and  off  (see  strong  oscillations  in  blue 
curve  in  Fig.  6). 

We  seek  a  non-dimensional  scaling  for  the  dissipation  rate  e  in  the  surfzone.  Following  Terray  et 
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Figure  5:  Estimated  surfzone  cross-shore  diffusivity  kxx  (with  error-bars)  versus  (a)  772g  T~l  and 
(b)  Vrot  Lx.  The  fit  slopes  are  11.7  and  0.2,  and  r2  are  0.32  and  0.59  for  (a)  //2„  T~  1  and  (b) 
Vrot  Lx,  respectively. 


al.  (1996),  who  developed  a  scaling  for  open-ocean  whitecapping  breaking  waves,  the  dissipation 
is  non-dimensionalized  as 


Hsige 

dF/dx 


(3) 


where  dF/dx  =  d{Ecg)/dx  is  the  cross-shore  gradient  of  the  the  incoming  wave  energy  flux 
(Ecg)  which  is  estimated  from  the  frames.  The  non-dimensional  e  has  a  consistent  relationship 
with  the  non-dimensional  distance  below  the  mean  surface  z' /Hsi g  at  each  frame  (Fig.  7).  This 
relationship  can  be  modeled  as  a  power  law  relationship  as 


Hsig£ 

dF/dx 


(4) 


where  C  and  7  are  fit  at  each  frame. 

At  the  3  surfzone  frames,  the  fit  7  vary  between  -1.5  and  -2  with  fit  skill  varying  between  0.4  and 
0.6  (Fig.  7).  This  best-fit  7  is  consistent  with  results  from  open-ocean  wave  breaking  which  has 
7  =  —1.9  (Terray  et  al.,  1996).  This  implies  that  the  mechanisms  by  which  turbulence  diffuses 
down  in  the  surfzone  are  similar  to  those  in  the  open-ocean.  However,  one  open  question  is  why 
the  constant  C  varies  so  much  between  the  frames.  This  is  being  investigated  further. 


IMPACT/APPLICATIONS 


Potential  impacts  include  improving  surfzone  and  nearshore  mixing  parameterizations  based 
upon  bulk  factors  such  as  wave  height,  wave  period,  bathymetry,  and  currents. 

RELATED  PROJECTS 
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Figure  6:  HB06  time  series  of  (top)  turbulent  dissipation  rate  e  at  frames  3  (blue)  and  6,  and 
(bottom)  incident  significant  wave  height  Hsig  at  frame  7. 


There  are  no  active  related  projects. 
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